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The measurement of elliptic flow for charged particles in Au-Au collisions at y/s NN = 130 
GeV using the PHOBOS detector at the Relativistic Heavy Ion Collider (RHIC) is presented. 
Charge particle distributions over a wide pseudo-rapidity region (\rj\ < 5.3) are measured by the 
PHOBOS silicon multiplicity detectors. These distributions are used to extract the final-state 
azimuthal anisotropy. The centrality and pseudo-rapidity dependences of elliptic flow, averaged 
over momenta on an event-by-event basis and particle species, are discussed. 



1. Introduction 

Collective flow in ultra-relativistic heavy ion collisions has become an important tool 
for searching for a hadronic phase transition since it is sensitive to the early stages of 
the collision system where the phase transition might occur, and its magnitude provides 
information on the Equation of State and the degree of thermalization of the system. 
In addition, the reaction plane information, a by-product of the flow measurement, will 
allow a better understanding in other physics measurements, such as HBT correlation or 
particle spectra. 

Flow is quantified as the coefficients of the Fourier expansion of the azimuthal angle 
distribution of particles with respect to the reaction plane: dN / d{4> — ^ r) = iVo(l + 
J2 n 2v n cos[n(</> — ^ /?)]), where denotes the reaction plane angle, and the first two 
coefficients, v\ and t>2, are the strengths of the directed flow and the elliptic flow, respec- 
tively. Both types of flow have been observed and extensively studied in lower energy 
experiments 
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2. Detector and data set 

Details of the PHOBOS detector can be found in Ref. [[|. The parts of the detector 
relevant for this analysis are, two sets of sixteen scintillator paddle counters (Paddle) for 
the event trigger and centrality measurement, two silicon tracking devices (Spectrometers) 
used for the event vertex reconstruction, the octagonal silicon pad detector (Octagon) and 
the annular silicon pad detectors (Rings) used to determine charged particle multiplicity. 

Collision events were selected by a trigger based on the coincidence of two or more hits 
in both Paddle counters. The signal in the Paddle counters increases with centrality and 
therefore can be used as a measure of event centrality. Further details on the trigger and 
on the centrality determination can be found in Ref. 

Reconstruction of the event vertex, i.e. the beam-beam collision point, starts with 
straight-line tracking based on a simple \ 2 fitting method in the first six planes of both 
Spectrometers where the magnetic field is negligible. The event vertex is then determined 
as the common vertex of the straight line tracks. 

Once the event vertex is known, hit distributions in the Octagon and Ring detectors 
are translated into particle pseudo-rapidity and azimuthal angle distributions. For colli- 
sions near the nominal interaction point, the Octagon and Rings cover |?7| < 5.3. Hits 
are defined as pads having more than 0.75 MIP[] energy after corrections for the parti- 
cle's incident angle and for the pad thickness. The definition of a hit is important as 
it changes the particle distribution itself in terms of both the total number of particles 
and the secondary/background ratio. Since a single particle can produce hits in several 
pads, especially in high rj region in the Octagon, an rj-<j) cell map of size 0.05x0.196 was 
introduced to allow merging of hits in that region. Finally, an occupied cell in the r]-(j) 
map is considered as a hit for this analysis. 



The analysis presented here is based on the data taken during the first RHIC physics 
run in 2000. 1.3 million events were accepted after checks of the detector performance 
and the vertex quality cut. A narrow fiducial cut was applied on the Z vertex position 
(-38cm< z vtx <-30cm) in order to ensure that the Octagon covers —2 < rj < 2 with a 
uniform and symmetric acceptance. The final sample used in the analysis amounts to 
13K events. The selected data are divided into 7 bins of centrality using the paddle mean 
signal as shown in Fig. [I]. 




Figure 1. The distribution of the mean 
Paddle signal normalized by the maximum 
mean Paddle signal for triggered events 
(upper curve) and for data accepted in the 
final analysis (lower curve). The dashed 
lines delineate the centrality bins used in 
the analysis. The numbers on the lines in 
each bin represent the integrated fraction 
of triggered events. 
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3. Analysis method 

The formalism for the flow analysis presented here is described in detail in |4[]. First, 
the event plane is reconstructed by the anisotropic flow itself by measuring the final state 
flow vector, Q = (X n ,Y n ) = w % cos(n0j), J2i w i s ^ n ( n 4>i))y where the sum runs over all 
particles, i, in an event and Wi are the weights to optimize the event plane measurement. 
The event plane angle \l/ n is then simply given as tan~ 1 (X n /y n )/n. Second, the observed 
flow signal is calculated as the coefficient of the Fourier expansion with respect to the 
event plane, v^ s =< cos(n(0 — \I/ n )) >. The actual flow signal v n measured with respect 
to the real reaction plane is then evaluated by taking into account the resolution in the 
event plane determination: v n = v°^ s / < cos(n(\I/ n — ^r)) >■ The event plane resolution 
can be directly calculated from the data using a sub-event technique since the correlation 
between two equal multiplicity sub-events yields a relation of < cos^^ 6 — ^r)) >= 
J< cos(n(^^ — ^n)) where a and b denote each sub-event. 

In this analysis, sub-events are defined by dividing particles into two pseudo-rapidity 
regions: one from -2.0 to -0.1 and the other from 0.1 to 2.0. The pseudo-rapidity gap 
between the two sub-events is introduced to reduce the short-range correlations between 
the sub-events. Detector related biases such as cracks between sensors, variation in 
acceptance of pads, and asymmetric production of background particles, are removed by 
taking the pad weights, w, as the inverse of the average particle density at each pad 
integrated over all events. After applying the weight correction, the \I/ n distributions 
become uniform. Fig. ^] shows the correlations between the sub-event plane angles. It can 
be seen that this analysis is quite sensitive to an elliptic flow v<i measurement, while the 
statistics of the sample are insufficient for a meaningful determination of v\. 
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Figure 2. Correlations between two event 
plane angles determined from the two sub- 
events: one from the negative pseudo- 
rapidity sub-event and one from the posi- 
tive pseudo-rapidity sub-event. 
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The second Fourier coefficients v™ w ( 



-< cos(2(0 a — ^/ b )) >) are evaluated by the az- 
imuthal angle distribution of hits in one sub-event (0 a ) with respect to the event plane 
from the other sub-event We use the notion of f™ 1 " because it is extracted from 

the 7]-(j) cell map, not directly from the particle angular distribution, and therefore this is 
not necessarily equal to v^ s as defined above. Monte Carlo (MC) simulation studies were 
performed to understand how v™™ is correlated to v^ s . It was found that v™ w decreases 
relative to v^ s due to the saturation of the detector with hits. The degree of the t>2 signal 
suppression, v^/vy*, could be sensitive to the detector occupancy and the magnitude 
of the flow itself. Extensive studies of simulated data show that the suppression has a 
negligible dependence on the size of flow for t> 2 < 10%, and can be parameterized by 
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1 — Occ), where Occ is the fraction of occupied cells in the rj-(j) map. Thus the observed 

u™/(l — Occ). The final elliptic flow strength, t> 2 , is 



flow value is deduced as v^ s 



derived by correcting v^ s with the sub-event plane resolution as described previously. 
4. Results 
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Figure 3. Elliptic flow, t>2, in the region 
-1.0< T] <1.0 function of centrality. 
The error bars are statistical only. 
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Figure 4. Elliptic flow, averaged over cen- 
trality, cLS cL function of rj. The error bars 
are statistical only 



The centrality and pseudo-rapidity dependences of v 2 are shown in Fig. |3] and Fig. £|, 
respectively. A strong centrality dependence of t>2 ranging from 0.03 for the most central 
events up to 0.07 for peripheral events is observed. The centrality dependence generally 
follows the results of a Hydrodynamic Model calculation. This is a possible indication 
that a high degree of thermalization is achieved in the early collision system by a large 
number of rescatterings. The present experimental result is in excellent agreement with 
an earlier measurement from the STAR Collaboration . The measurement of v 2 at high 
pseudo-rapidity, showing a similar shape to dN/dr], is a unique result from PHOBOS. 

Various cross-checks were done using MC simulation to confirm the analysis method by 
studying the consistency between input flow and the resulting value. An overall systematic 
error of 0.007 is estimated from studies of numerous sources of systematic uncertainty: (1 — 
Occ) parameterization, hit definition, beam orbit, relative alignment between Octagon, 
Rings, and Spectrometers. Systematic errors related to the vertex algorithm, weighting 
matrix, sub-event technique, and magnetic field settings, were all found to be negligible. 
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